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ABSTRACT 

We present a near-infrared spectroscopic study of a stellar mass selected sample 
of galaxies at z ~ 1 utilising the LIRIS multi-object spectrograph on the William 
Herschel Telescope. We detect continuum, and the Ha line for our sample, which is one 
of the better direct tracers of star formation in external galaxies. We spectroscopically 
measure the Ha emission from 41 massive (M, > 10^" '^ M© ) galaxies taken from 
the POWIR Survey with spectroscopic redshifts 0.4 < Zgpcc < 1-4. We correct our 
Ha fluxes for dust extinction by using multi-wavelength data, and investigate SFR 
trends with mass and colour. We find a drop in the fraction of massive galaxies with 
M* > 10"'^^ Mq which are detected in Ha emission at z < 0.9. We furthermore find that 
the fraction of galaxies with Ha emission drops steadily and significantly with redder 
{U — B) colours at z ^ 1, and that the SSFR drops with increasing [U — B) colour for 
galaxies at all masses. By investigating the SFR-mass relation we find that the SFR is 
roughly constant with mass, in possible contrast to previous work, and that the specific 
star formation rate (SSFR) is lower in the most massive galaxies. The scatter in the 
SFR vs. mass relationship is very small for those systems with ongoing star formation 
which suggests that star formation in the most massive galaxies at z ~ 1 shuts off 
rather abruptly over < 1 Gyr, without an obvious gradual decline. We furthermore 
investigate the SFR as a function of (C/ — B) colour divided into different mass bins, 
revealing a tracer of the epoch of transition from star forming to passive, as a form of 
star formation "downsizing" . This suggests that the shut off of star formation occurs 
before the change in a galaxy's colour. We find that galaxy stellar mass is the primary 
driving mechanisms behind the star formation history for these galaxies and discuss 
several possible mechanisms for regulating this process. 

Key words: galaxies: star formation - galaxies: evolution - galaxies: formation 



1 INTRODUCTION 

It has been long known that galaxies have evolved over the 
history of the universe, and a major branch of recent as- 
tronomy has been to quantify how the galaxies we see in 
the present day came to be as they are now. It is impossible 
to observe how a given galaxy evolves over cosmologically- 
significant time scales, but its evolution can be inferred ei- 
ther on a galaxy by galaxy basis from the present state of 
individual galaxies - for example the star formation his- 
tory (e.g. Maraston et al. 2006; Noeske et al. 2007b), or the 
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merger history (e.g. Conselice 2003; Conselice et al. 2009); 
or from the changes to the properties of a population of 
galaxies at different times. 

Massive galaxies are a major component of a number 
of branches of astronomy, being a pivotal test for, amongst 
others, dark matter theories and galaxy formation models. 
Multiple constraints have now been put on when the most 
massive (M« > 10^" Mq ) galaxies formed (e.g., Fontana et 
al. 2004; Glazebrook et al. 2004; Bundy et al. 2006; Conselice 
et al. 2007, 2011) and it is seen that 2; ~ 1 is an era of tran- 
sition when massive galaxies shut off their star formation; 
their mass for the most part, already in place (Glazebrook 
et al. 2004; Bundy et al. 2006; Conselice et al. 2007; Perez- 
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Gonzalez et al. 2008a; Ilbert et al. 2010). The star formation 
in these most massive galaxies was significantly higher be- 
fore this epoch, and there appears to be evidence of star 
formation rates well in excess of those found in today's uni- 
verse (e.g. Hopkins 2004; Perez- Gonzalez et al. 2008b). 

Within this transition epoch, there is increcising evi- 
dence that the shut off in star formation is not triggered si- 
multaneously across the universe, but is highly dependent on 
the masses of the galaxies. In a process that has been named 
'downsizing', it has been suggested that the most massive 
galaxies finish their star formation earlier than less massive 
systems, a progression observed in broad-band colours (e.g. 
Bundy ct al. 2006), radio observations (Hopkins 2004), in- 
frared observations (e.g. Papovich et al. 2006; Bell et al. 
2007) and emission line surveys (e.g. Heavens et al. 2004; 
Panter et al. 2007). 

Measurements of high-redshift star formation rates 
(hereafter SFRs) have until recently, been dogged by many 
problems. The more favoured emission line diagnostic, the 
Ha line is only available via optical spectroscopy at z < 
0.4 and so the more troublesome [Oil] line has been regu- 
larly used at higher redshifts (Gallagher, Hunter & Bushouse 
1989; Kennicutt 1998a; Rosa-Gonzalez, Terlcvich & Tcr- 
levich 2002; Chariot et al. 2002; Kewley, Geller & Jansen 
2004). The [OH] line is affected more by metalicity and due 
to its bluer wavelength, is more influenced by dust atten- 
uation than the Ha line (Jansen, Pranx & Fabricant 2001; 
Chariot et al. 2002). As a result, it is never entirely certain 
whether high-2 SFRs derived from [OH] are completely com- 
patible with \ow-z SFRs derived from other more accurate 
methods. 

Infrared (IR) SFR tracers, for example 24 /im emission 
(Lonsdale Persson & Helou 1987; Rowan-Robinson & Craw- 
ford 1989; Kennicutt 1998a; Calzetti et al. 2007) is also com- 
monly used as a 'proxy' SFR diagnostic. Emission from SFR 
regions is absorbed by dust and re-radiated as IR radiation, 
allowing a conversion to be calibrated to retrieve the origi- 
nal SFR. This method has however, problems in measuring 
star formation rates. For example, using current technol- 
ogy, spectral line tracers are also able to measure SFRs to 
an order of magnitude lower than the Spitzer MIPS IR in- 
strument, and are more sensitive at higher redshifts (see for 
example Noeske et al. 2007a). Other problems include the 
fact that the extrapolation from a few sub-mm and radio lu- 
minosities to a bolometric IR luminosity is subject to many 
uncertainties, especially at high rcdshift (Small, Ivison & 
Blain 1997; Chapman et al. 2005) and systematic changes 
(for example Papovich et al. 2007); contamination by an 
active galactic nuclei (AGN) (sec Daddi ct al. 2007b); and 
possibly cirrus emission comming from older stars (Lonsdale 
Persson & Helou 1987; Helou 1986; Kennicutt 1998a). Radio 
SFRs, another SFR measure, often relies on stacking of data 
sources which may lead to complications in interpretting re- 
sults (e.g. Ivison et al. 2010). 

Each measure of star formation tracks a different cle- 
ment of the stellar evolution cycle, each with its own precise 
definitions (See Kennicutt 1998a for a review of what dif- 
ferent wavelengths actually measure) and therefore unique 
biases. It is important therefore to get a measure of whether 
the apparent changes in global SFR are actually due to a 
change in the SFR of galaxies, or because the different mea^- 
surements at different redshifts are measuring different pro- 



cesses. The Ha line is the best tracer of this due to the 
relatively low attenuation it suffers from dust compared to 
other emission line tracers, the lack of any strong depen- 
dency on metalicity, and because its emission is associated 
largely with only star forming regions. It is also a direct mea- 
surement of star formation, and measures instantaneous star 
formation during the last 20 Myr within the galaxy. Atten- 
uation is still a significant problem however - it is variable 
galaxy-to-galaxy due to its dependence on the amount of 
dust in a galaxy, and can effect the emission by a factor of 
three or more. 

At redshift z > 0.4, the Ha emission line shifts into 
the near-Infrared (NIR) and different technologies arc re- 
quired to observe it. In the last decade there have been 
a rmmber of IR instruments built for telescopes and with 
these new instruments, we can now begin to measure this 
emission to higher redshifts than before. There have been 
a few small investigations into Ha emission at z ~ 1, how- 
ever most have been hampered by a very high detection 
limit, or the low number of galaxies studied. Glazebrook 
et al (1999) presented the first NIR spectroscopy to detect 
Ha emission from galaxies at 2 = 1. They studied 13 field 
galaxies drawn from the Canada-France Redshift Survey and 
compare their derived SFRs with observations of the rest- 
frame UV emission. They find that the Ha SFR density 
measurements are roughly three times larger than the ones 
derived using the Madau, Pozzetti & Dickinson (1996) UV 
relation. Dohcrty et al (2006) followed from works hy Yan 
et al. (1999), Hopkins, Connolly & Salay (2000) and Tresse 
et al. (2002) examining the SFR density at 2 ~ 1. Doherty 
et al. took spectra of 38 galaxies within the Hubble Deep 
Field North and stacked their spectra in order to measure 
the total star formation and hence the SFR density. They 
estimate a lower limit for the total SFR in their stacked 
sample to be ~ 312M0 yr~^ but do not measure individ- 
ual SFRs. Both Doherty et al. and Yan et al. derive dust 
corrections by comparison with UV measurements. Tresse 
ct al. (2002) used a mixture of Balmer ratios (Ha/H/J and 
H7/H/J) where available to measure a dust correction, in 
other galaxies they used a fixed attenuation. Most recently, 
Rodn'gucz-Eugonio ct al. (2007) measured 16 star forming 
galaxies at redshift 0.8 < z < 1.0 from the DEEP2 survey 
with the LIRIS camera on the William Herschel Telescope, 
but had problems removing the sky lines from their data. 
Garn et al. (2010) used narrow band imaging to measure 
Ha emission from 477 sources in the HiZELS survey with a 
SFR detection limit of 2.15 M© yr~^. They used other avail- 
able multiwavelength imaging to calculate the dust attenua- 
tion using a number of methods. Sobral et al. (2010) further 
investigate the SFR-mass relation and specific SFR-mass re- 
lation of these galaxies, however focusing on environmental 
effects. 

None of these previously mentioned surveys have how- 
ever used Ha spectroscopy to investigate trends and rela- 
tions found with other SFR tracers, nor do they use a mass 
selected sample. We present the first such study using the 
LIRIS spectrograph on the William Herschel Telescope to 
obtain near-Infrared multi-object spectroscopy of 42 galax- 
ies from the POWIR survey (Consclice et al. 2008a) with 
stellar masses M* > 10^" and at redshifts 0.4 < z < 1.3. 
We measure Ha luminosities for our sample to calculate 
their SFRs and to investigate whether the observed trends 
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in global and individual star formation properties are seen 
when using different SFR tracers. In this paper we further- 
more investigate the number of galaxies with Ha emission 
as a function of mass, compare Ifa-derived SFRs with other 
measures of their SFRs and relations between the SFRs, and 
mass and colour. 

This paper is organised as follows. The data and new 
observations are summarised in Section[2l including different 
methods of correcting for dust attenuation in Section 12.51 
The galaxies where Ha emission is detected are investigated 
in Section|3]and the SFR relations are investigated in Section 
13.51 The results, and the astrophysical interpretation of the 
SFRs are discussed in Section[4l with a summary of the work 
in SectionO Throughout the paper, we adopt the cosmology 

= 0.3, VIa = 0.7 and Ho = 70kms"^Mpc"\ 
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2 DATA & SAMPLE 
2.1 Sample 

The galaxies investigated in this paper are selected from 
galaxies within the POWIR survey (Conselice et al. 2008a) 
as investigated in Conselice et al. (2007, hereafter C07). The 
POWIR survey is a large area (~1.5deg^) deep NIR survey 
in the K and J band covering the GOODS field North (Gi- 
avalisco et al. 2004), the Extended Groth Strip (EGS; Davis 
et al. 2007) and three fields observed by the DEEP2 team 
with the DEIMOS spectrograph (Davies et al. 2003). C07 
used all the areas covered except the GOODS-North field. 

Data in C07 was acquired from the Palomar 5-m tele- 
scope [Ks and J band imaging) and used imaging from 
the Canada-France-Hawaii Telescope (CFHT; B, R, I), 
the Advanced Camera for Surveys (ACS) on the Hubble 
Space Telescope (Optical), the Multiband Imaging Photome- 
ter (MIPS) on the Spitzer Space Telescope (24 /im), the 
Very Large Array (1.4 GHz), the Chandra X-ray Obser- 
vatory and spectra from the Deep Imaging Multi-Object 
Spectrograph (DEIMOS) at the W. M. Keck Observatory 
(A = 6500 - 9100A). A .K'-selected catalogue was created 
and the different photometry bands matched. Readers are 
directed to C07 for more information on these observations 
and the parent sample selection. 

The C07 catalogue was used to measure the stellar 
masses of galaxies to if ~ 20 - 21 from SED fitting as 
in Bundy et al. (2006; for details, see G07). This measure- 
ment uses a Chabrier initial mass function (hereafter IMF; 
Chabrier 2003). A complete subset was then created with 
M* > 10^° Mq and a redshift range of z ~ 0.4 - 1.4 taken 
from DEEP2 spectroscopic redshifts (Davis et al. 2003). This 
is the dataset from which our sample was constructed. 

The targets for LIRIS were taken from two regions, 
one purely mass selected from the POWIR Field 3 to have 
M. > 10^° Mq , the second selected from the EGS field. This 
second set of targets were selected to be massive (M, > 
10^° Mq), have a significant MIPS flux and no Chandra 2-8 
keV X-ray emission at < 10*^ ergs~^ cm~^. Both sets were 
selected within the range 0.4 < Zspcc < 1.3, giving a total of 
~7500 possible targets. Masks for multi-object spectroscopy 
were created using pointings that maximised the number of 
galaxies targeted. 



Figure 1. The redshifts and stellar masses of the POWIR survey 
(grey points) and our observed galaxies taken from this sample 
(black triangles). The dashed line is our mass cut. 



2.2 NIR Spectra 

Our Near infrared (NIR) spectra were taken during 9 nights 
over two observing runs on the 4.2-metre William Herschel 
Telescope using the Long-slit Intermediate Resolution In- 
frared Spectrograph (LIRIS) instrument in Multi-Object 
Spectroscopy mode (Acosta-Pulido et al. 2003). The spec- 
tra were taken using a 3-point nodding routine and the total 
integration times range from 2.1 hrs to 6.3 hrs for each of 
the 7 masks. Each mask targeted between 8 and 14 galaxies, 
targetting a total of 71 galaxies. 

Data reduction was done using the dedicated suite of 
routines lirisdr for iraf, written by Jose Acosta-Pulido 
and the spectra were flux calibrated using J-band magni- 
tudes from the POWIR survey data. In all cases the spectra 
covered the wavelength range of the J-band, and the spec- 
tra were scaled to give the same integrated flux measurement 
within that wavelength range as the J-band magnitude mea- 
surements. This also corrected for flux lost due to the flnite 
width of the slit (aperture correction). Out of the 71 galaxies 
targeted, 42 had detected continuum. There are no galaxies 
where there is emission at around 6563 A but have no de- 
tected continuum. There were 21 galaxies detected in each 
of the two areas of sky targeted. One galaxy, 32014349 is af- 
fected by an atmospheric absorption feature within 200A of 
the wavelength of any Ha emission line and so is discarded 
leaving 41 galaxies in our sample. The stellar mass and red- 
shifts of these 41 galaxies are shown in Figure [1] (see also 
Table©. 

The ID Spectra were examined by eye and if an emis- 
sion line was seen at the correct wavelength to be Ha, the 
flux was extracted from the continuum subtracted spectra 
within this region. Even though in most cases where there 
was substantial Ha emission the HaA6563 and [NIIJA6583 
appeared separated, the extraction region was chosen to in- 
clude both emission lines. This was to allow the procedure 
for measuring an upper limit where there was no emission 
detected to follow as closely as possible the procedure used 
when there was emission. The continuum level and noise was 
taken from two 200 A sections of continuum one average full 
width half maximum (FWHM, 5 pixels = 31 A) either side 
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of the emission and the root mean squared value of this 
continuum is used as the error on each pixeL This average 
FWHM is used as the extraction width for continuum with 
no visual signs of emission. 

If the Ha emission is less than the 1.5 times the er- 
ror and no line has been seen by eye, then these galaxies 
are classified as non-detections. This method of measuring 
the emission line is similar to that used in Glazebrook et al 
(1999). A correction of 30 per cent for [Nil] is applied (see 
for example, Gallego et al. 1997; Tresse et al. 1999; Garn et 
al. 2010). In cases where there was no visible emission, the 
continuum noise was used to calculate a 1.5 a upper limit to 
any Ha emission. The DEEP2 spectroscopic redshifts were 
used to determine the predicted observed Ha wavelengths 
and convert fluxes to luminosities. 

2.3 Star formation rates 

Star formation rates were derived from the Ha luminosi- 
ties using the Kennicutt (1998) relation converted to use a 
Chabrier IMF (Chabrier 2003): 

SFR (M0 yr"') = 4.6 x IQ'*'^ L{Ha.) (ergs"'). (1) 

In order to compare our measurements with measurements 
taken in other wavelengths, we use Infrared (IR), UV and 
radio observations already available of the same sources. Ul- 
traviolet (UV) data was taken using the Galaxy Evolution 
Explorer (GALEX) and available imaging in the B, R and / 
filters to match the galaxy's rest-frame emission as close as 
possible to the GALEX filters (1528 A and 2271 A; Schimi- 
novich 2007). SFRs and dust extinction-corrected SFRs were 
calculated using a composite method using IR calibrated 
UV-optical colours and broadband SED fitting (for a de- 
scription of the method see Schiminovich et al. 2007 and 
Seibert et al. 2005). The UV and IR SFRs are discussed in 
more detail in Section [2.51 

Radio data was acquired from a 1.4 GHz catalogue of 
the EGS, (Ivison et al. 2007) taken with the Very Large 
Array (VLA). They detect 1,123 sources with S^^qjj^ 
50 /xJy to a 5-a detection limit which corresponds to rest- 
frame 1.4 GHz calculated SFR of 275 Mq yr'^ at z = 1.0. 

Our sample were matched with the radio catalogue, but 
only reveal 1 match (from 41 galaxies). The radio catalogue 
was then matched with the entire set of possible galaxies 
that could have been observed and are in the same region 
of sky as those that were observed. Out of these 62 galaxies 
there are 4 matches (6 ± 3 per cent). 

Concluding that our galaxies are not detected at 1.4 
GHz, upper limits can be placed on the total star formation 
in these galaxies following the methods of Condon (1992), 
Haarsma et al. (2000) and Condon, Cotton & Broderic 
(2002). Using the detection limit of the nearest galaxy in the 
sample (z ~ 0.55) gives an upper limit of 66 yr~^, the 
furthest (z ~ 1.23) gives ~450Mo yr"\ These are consid- 
erably higher than the majority of our upper limits derived 
from the Ha measurements and thus are not further used. 

2.4 Imaging 

Figure [2] shows the subset of our sample imaged by the Hub- 
ble Space Telescope ACS camera through the F606W filter 



with the position of the LIRIS slits overlaid. In general, the 
slits are well aligned with the target galaxies. For all galax- 
ies, the slit captures the centre and some of the outer-regions 
and so we can be confident that we are capturing a reason- 
able representation of the emission from star formation. 

The majority of the galaxies in our sample are clearly 
disk galaxies. There are a couple of galaxies with disturbed 
morphologies and signs of possibly being merging systems. 
Note that none of the galaxies in our sample were selected by 
morphology. Below we list some features of our four systems 
whose observations may be contaminated by nearby galaxies 
or stars. 

Galaxy 12020027 (Figure [2};) is actually two galaxies 
separated by 1.4". These two galaxies are close enough to- 
gether that the 24-/im infrared data from the Spitzer Space 
Telescope will not differentiate between the two galaxies due 
to its large point spread function (PSF). If the second galaxy 
has a higher dust content than the target galaxy, the Ha 
emission will be dominated by the galaxy with the least 
dust, whereas the 24-/im emission is dominated by emission 
from the dusty galaxy. This will lead to an applied dust 
correction that is too large (or a too large amount of Ha 
emission being corrected). This is a possible explanation for 
the high dust-corrected SFR from this galaxy (see Section 
[23]) . 

Galaxy 12020031 (Figure [2ji) has a possible second 
galaxy 2.0" away from the centre of the main galaxy. Again 
this may affect the dust-corrected SFR measured. 

Galaxy 12024436 (Figure (2)) appears to have a small 
object, probably a foreground star also in the slit, 1.5" away 
from the galaxy. This is not an obvious extra component 
in the spectra, but may affect other measurements of the 
galaxy. This galaxy has the appearance of a merging system 
with two central peaks separated by 0.6" 

Galaxy 13026831 (Figure [2J1) also appears to have two 
parts separated again by 0.6" 

2.5 Dust corrections of star formation rates 

Each method of measuring the star formation rate (SFR) of 
galaxies uses emission at certain wavelengths. Optical trac- 
ers of SFR are affected by dust extinction that is a function 
of wavelength, so we must apply these corrections. In the 
NIR, optical, and UV wavelengths, the effect is manifested 
as an attenuation to the emission that is highly dependent 
on wavelength (see for example, Calzetti 2001). 

Many studies - usually at redshifts where observations 
of emission lines becomes difficult, or where spectroscopic 
data is not yet available - often use a blanket value for the 
amount of attenuation a galaxy suffers at a specific wave- 
length, A{\) (often 1 magnitude at V-band: A{Ha) ~ 0.8) 
and apply it to all the galaxies in their sample (e.g. Kenni- 
cutt 1983; Chariot et al. 2002; Tresse et al. 2002; Garn et al. 
2010 and references within). There is evidence however, that 
there are a range of attenuations in galaxies at 2; ~ 1 (e.g. 
Garn et al. 2010). In addition, there is also a large uncer- 
tainty of up to ~30 per cent on the conversion of Ha lumi- 
nosity to SFR due to different values derived from different 
stellar population models (see Kennicutt 1998a). Therefore 
this method can crudely work to statistically correct a large 
sample. Our sample of galaxies is neither large enough, nor 
statistically general enough to apply such a method. 
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a) 12019914 b) 12019916 c) 12020027 d) 12020031 




i) 12024427 j) 12024436 k) 12024440 1) 12024445 




Figure 2. Hubble Space Telescope ACS F814W images of the galaxies where HST data is available in the EGS field (Conselice et al. 
2007) overlaid with the slit position and width (white lines). The field of view of ea<;h image 9" by 9". 
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0.79 


1.01 


10.58 


<0.9 


0.60 


<0.7 


<3.0 






32008Q1 5 


0.74 


0.75 


10.56 


2.1 it 0.9 


0.61 


1.6 ± 0.7 


6.7 ± 9.7 






32008Q24 


0.78 


0.83 


10.96 


6.5 it 1.5 


0.63 


5.1 ± 1.2 


22 1 + 27 2 






3200Q330 


0.90 


1.05 


11.01 


<4.2 


0.62 


<3.3 


<13.9 






32014288 


0.78 


1.03 


10.83 


<3.9 


0.61 


<3.0 


<12.8 


- 


- 


32014399 


0.79 


0.66 


10.63 


50.9 it 3.6 


0.81 


40.2 ± 2.9 


266.7 ± 285.7 






32014580 


0.79 


0.61 


10.46 


15.6 it 2.4 


0.66 


12.3 ± 1.9 


57.9 ± 66.9 






32018834 


0.81 


0.59 


10.08 


2.9 ± 1.3 


0.61 


2.3 ± 1.0 


9.5 ± 13.7 






32018902 


0.82 


1.29 


10.72 


<7.3 


0.63 


<5.8 


<25.0 






32018960 


0.78 


1.14 


11.03 


<18.5 


0.68 


<14.6 


<70.7 






32019661 


0.78 


1.02 


11.03 


<3.0 


0.61 


<2.3 


<9.7 






32019921 


0.78 


1.17 


10.99 


<5.8 


0.62 


<4.6 


<19.6 






32019951 


0.80 


1.05 


11.02 


<9.7 


0.64 


<7.6 


<33.8 






32024002 


0.74 


1.31 


10.71 


<11.5 


0.65 


<9.1 


<41.2 






32024419 


0.79 


0.65 


10.43 


<11.3 


0.64 


<8.9 


<40.2 







Table 1. Galaxies observed in our sample. (U — B) is colour, logM, is the stellar mass, Ljja is the Ha luminosity, E{B — V) is the (B — V) 
colour excess, SFRHa,uncorr! SFRjja,corr) SFR24 /im and SFRuv are the star formation rates measured via Ha emission uncorrected and 
corrected for dust, 24-;um emission and UV emission. The Ha derived SFR is not corrected for dust, the UV derived SFR has been 
corrected. 



There are two sources of emission used at these wave- 
lengths for measuring the SFR nebula emission lines (e.g. 
Hq, H/3, Paa, Oil) and stellar emission (usually in the UV). 
Emission lines are caused by recombination or coUisional ex- 
citation in HII regions which are formed by massive, short- 
lived O and B type stars; similarly the UV emission gives 
the equivalent of a number count of the these stars. Both of 
these have been calibrated with certain assumptions to give 
simple conversions to the SFR of the galaxy (see Kennicutt 
1998a). 

The attenuation of a galgixy at a specific wavelength, 
A(X), is related to the colour excess of the nebular emission 
in the galaxy, E{B — V)gas as given by: 



A{X) = E{B - y)ga. fe(A) (2) 

where k{\) is the value of the dust extinction curve be- 
ing used at the wavelength of emission. There is signifi- 
cant large-scale differences between different dust extinction 
curves that have been measured to date. However a study by 
Calzetti (2001) showed that between the (rest-frame) wave- 
lengths of ~3700 A and ~6600 A the differences are minor. 
The Galactic extinction curve of Cardelli, Clayton & Mathis 
(1989) is used in this study which gives k{Ha) = 2.54, 
k{Hl3) = 3.61 and k{H-f) = 4.17. 

A standard method to correct for dust on a galaxy by 
galaxy basis is to measure the ratio of the Ha to H/3 emis- 
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ID 


EiB-V), 


;as from 




H«/H/3 


H7/H/3 


12019914 




3.8 ± 1.3 


12024309 


1.98 ± 0.43 


-0.54 ± 0.64 


12024427 




2.60 ± 0.72 


12024440 


0.83 ± 0.69 


0.83 ± 0.99 


12024445 


3.8 ± 1.6 




12024528 


0.91 ± 0.49 




13027461 




0.30 ± 20 


13035323 


-0.08 ± 0.50 


0.57 ± 0.22 


13035770 




1.59 ± 0.60 


13035775 




2.1 ± 1.1 



Table 2. E{B — V)ga,B values derived from the Balmer decrements 
for the galaxies where line fluxes are available. Values would be 
expected to be — V)gas ^ 1-5; negative values are unphysical 



sion and compare this to the theoretical value (the Balmer 
decrement). As these lines will be attenuated by different 
amounts, a dust extinction curve can be used to calculate 
the amount of attenuation. This method can be used with 
any two emission lines that have a known intrinsic ratio. Five 
galaxies in our sample have both Hq measurements (not 
upper limits) from our measurements, and H/? luminosities 
from the DEEP2 spectroscopy allowing for a calculation of 
the dust correction through the Balmer decrement. Eight 
galaxies have H7 from DEEP2 allowing a comparison with 
the H/3 lines. Three galaxies have all three values. Intrinsic 
ratios of {Fi^'/Ff'^) = 2.87 and {Ff'^/Ff'^) = 0.466 (Case 
B recombination at T = 10,000K, Osterbrock 1989) are used 
to calculate E{B — V)gas values. 

Even within the large errors, we find significant differ- 
ences in the range of measured values (Table [2]), so that un- 
fortunately this extinction correction method does not give 
satisfying results. The large errors are likely the result of in- 
sufficient signal-to-noise ratio of both our present Ha survey 
and the DEEP2 survey's H/3 and H7 measurements. 

We can also use emission in the far infrared to calculate 
dust corrections, as UV light absorbed by dust is re-radiated 
in the far-IR and this can be used to correct the attenu- 
ated SFRs measured in other wavelengths. Data taken from 
the Spitzer MIPS camera was if-corrected and converted to 
SFRs by Noeske et al. (2007). As this is a measure of ob- 
scured SFR, this data can then be combined with UV data 
to correct for dust. 

We use two methods to account for dust. A direct but 
simplistic method to obtain the total SFR is to use a summa- 
tion of the un-obscured SFR from the UV and the obscured 
SFR from the 24^m. 

SFRxot = SFRuv + SFR24 Mm , (3) 

where, when calculating the SFRuv from the far-UV (FUV) 
flux (flux detected at the wavelength of the GALEX FUV 
band which has an effective wavelength of 1528 A) 

SFRuv (Mq yr"') = 10-'*-'''"Lpuv(ergs-i Hz"!), (4) 

(Salim et al. 2007) and the SFR calculated from the IR 
emission 



SFR24m™ (M0 yr-1) = 

1.27xl0-^« [L24Mm(ergs-i)]™, (5) 

(Calzetti et al. 2007). This requires accurate calibrations of 
both the derived SFRs, and ignores possible effects from dust 
geometry and complex radiative transfer processes. Kenni- 
cutt (2009) derive a set of coefficients that are multiplied 
by an IR SFR tracer luminosity to convert it into the lumi- 
nosity from an attenuated emission-line SFR tracer that has 
been absorbed by dust. In both of these methods the IR SFR 
(luminosity) is much larger than the Ha SFR (luminosity) 
and therefore the corrected SFR becomes approximately the 
IR SFR 

Schiminovich et al. (2007) developed a method for cor- 
recting the UV-derived SFRs for dust attenuation using 
a hybrid approach combining an IR-calibrated measure of 
the far-UV (FUV) attenuation {Amx) based on UV-optical 
colours and the 4000A break Z)„(4000) from Johnson et al. 
(2006, 2007) (accurate for galaxies with D„(4000) < 1.7) 
and an attenuation measure that is based on detailed 
model fits to the SDSS absorption-line spectrum and broad- 
band SED from Kauffmann et al. (2003a) (likely to be more 
accurate for galaxies with higher Z)„(4000)). The combined 
fit is used with the Calzetti et al. (2000) extinction curve to 
correct the UV luminosity, and then the relation from Salim 
et al. (2007, Equation|3| to derive the SFR. For more details 
and discussion, see Schiminovich et al. (2007). We find that 
13 out of our sample of 41 galaxies have SFRs that can be 
derived from GALEX UV data. 

Using either of these methods, we can calculate the to- 
tal amount of star formation in each galaxy. This can then 
be compared with the SFR measured using Ha emission to 
determine the amount of extinction in the Ha fiux. As the 
UV emission comes directly from stars within a galaxy, the 
ratio of detected UV emission to dust corrected UV emission 
can be used to measure the E{B — V) value for the stellar 
component of the galaxy. This E{B — l/)steiiar value must 
be converted to the E{B — V)gas value for use with neb- 
ular emission lines. The relationship has been empirically 
measured as 

E{B - VU.iu. = (0.44 ± 0.03) E{B - F)gas. (6) 

(Calzetti, Kinney & Storchi-Bergmann 1994; Calzetti et al. 
2000). 

There appears to be a correlation between a galaxy's 
dust content and its SFR (e.g. Hopkins et al. 2001; Sullivan 
et al. 2001; Bauer et al. 2011). This may arise from the fact 
that the SFR in a galaxy is proportional to the gas surface 
density (Kennicutt 1998b) and so galaxies with higher SFRs 
must have more gas and dust and hence more extinction. 
From this idea, many studies have tried to create either a 
theoretical or empirical law that relates the amount of star 
formation in a galaxy to the dust extinction. Three such 
relations are described by Afonso et al. (2003), Choi et al. 
(2006) and Garn et al. (2010). 

Figure |3] shows the comparison between E{B — F)gas 
correction values derived from the UV data (Equation |31 
with the corresponding E{B — V^)gas values derived from 
the total SFR (Equation [3)1 and with the E{B - l/)gas val- 
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Figure 3. a) Comparison of tlie UV-derived E{B — V)gas values {E{B — V)uv) witli tlie E{B — y)gas values derived from the total 
UV+IR SFR {E{B - V)Tot)- The dashed line is the 1:1 relationship, b) Comparison of the UV-derived E{B - y)gas values with the 
empirically derived values {E{B — V)Empirical) from relations in Garn et al. (2010, crosses), Choi et al. (2006, triangles) and Afonso et 
al. (2003, diamonds). The relations in Choi et al. and Afonso et al. have been converted from a Salpeter IMF to a Chabrier IMF by 
lowing the SFR of each point by a factor of 1.7. 
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Figure 4. Comparison of the E{B — y)gas values derived from 
UV measurements as a function of uncorrected SFR with the 
empirically derived relations shown. The solid line shows the best 
fit relation from the UV-derived values. The three curved lines 
are the empirical models, Afonso et al. 2003 (dot-dot-dash), Choi 
et al. 2006 (dot-dash), Garn et al. 2010 (dotted). 

ues from the empirical relations of Alfonso et al., Choi et al., 
and Garn et al. The E{B — V)gas values derived from the 
UV are consistent with the values derived from the empir- 
ical methods. Figure 3] shows the UV-derived E{B — V)gas 
as a function of observed observed SFRhq, along with the 
empirical relations. There is a slight increase in E{B — V)gas 
with higher SFR, and this can be fit with a straight line by 
a least-square method function as: 

E{B - V)vY = 

(8.96 ± 13.7) X 10"^ SFRna.uncorr + (0.663 ± 0.113) , (7) 

Although with significant scatter. Each method for deriving 
the E{B — V)gas value for each galaxy has its own biases, 
short comings and advantages, and as there are such differ- 
ences between the methods, they cannot be interchanged. 
The median E{B — V^)gas value for each galaxy is calcu- 



Figure 5. Comparisons of the E{B — V)gas values derived from 
the UV SFRs with the median value of all methods. Most UV- 
derived values lie near the median value, or are lower 

lated from the E{B — V) values derived from the different 
available methods. However because not all the methods can 
be used for all the galaxies, the median value of E{B — V) 
for each galaxy will not be comparable across our sample. 
Some galaxies do not have enough data for any corrections, 
and where galaxies only have enough data for one method, 
it then introduces that method's specific biases onto that 
galaxy only. Picking only one method ensures the method 
based biases for all galaxies are the same and facilitates some 
attempt at correcting for this. 

We thus use the E{B — V)gas value derived from the 
GALEX UV data to correct our Ha-derived SFRs which al- 
lows the most galaxies to be corrected with their own ob- 
servational data. The fitted relation in Equation [7] is used 
to correct galaxies for dust extinction with no UV derived 
SFRs, which equates to an attenuation that ranges from 
Ahc = 1.68 at SFRHcuncorr = 0.1 M© yr"^ to Ahc = 2.25 
at SFRh a ,uncorr 

= 25 Mq yr-^ Fi gure [S] shows the devia- 
tion of the UV-derived E{B — V')gas values from the median 
E{B — V)ga.s value for that galaxy calculated over all our 
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Figure 7. The redshifts of galaxies with M» > 10^^ divided 
into those with Ho detections and those without. All detections 
of star formation in these massive galaxies are at z > 0.93 

methods. The UV-derived E(B — V^)gas values are seen to 
be in reasonably good agreement with the median values. 



3 RESULTS 

3.1 Ha detections 

Out of the 41 galaxies with a continuum detection in the 
sample, 14 (34 ± 10 per cent) have significant detected Hq 
emission. The median detection limit for uncorrected star 
formation is 3.0 ± 0.5 M© yr~^. First, to investigate which 
type of galaxies are being detected, the fraction of detections 
are plotted as a function of redshift in Figure[6K. At redshifts 
higher than z — 0.7, there is no statistical change in the 
detection rate as a function of redshift, implying there is 
no strong bias against galaxies at either end of our redshift 
range in our sample. 

The errors plotted for each point are the combination 
of the error on the number of galaxies detected in Ha, ctjv^ 
and the total number of galaxies, o-jv^, assuming Poisson 
statistics for both, with a cross term that takes into account 
the fact that these two quantities are not independent. 

(Bevington 1969). cr%^f^^ is the covariance between Nd 
and Nt and for a positively-correlated pair of gaussian- 
distributed sets, can be approximated as 

f^N^NT — '^Nj^t^NT- (9) 

Figure [SJj plots the detection fraction as a function of the 
stellar mass of our galaxies. There is a clear, significant and 
sharp drop off in the fraction of galaxies with detected Ha 
emission with M, > 10^^ Mq . In the most massive bin there 
are a total of 2 galaxies with Ha detections (12 ±25 per cent) 
and 14 with no detections. The redshifts of these galaxies are 
shown in Figure[71 and the galaxies detected in Ha are all at 
around 2 ~ 1. The distribution in mass of the non-detections 
is similar to the whole sample. Although our sample is small. 
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Figure 9. The averaged luminosity spectra for all galaxies with 
no Ha detected. 



there are no galaxies with M, > 10^^ Mq that have Ha de- 
tected at a redshift below z — 0.93 (Figure [7]). 

Our Ha detection fractions are plotted as a function of 
{U — B) colour in Figure [EJ;. Bluer galaxies have a higher 
fraction with detected Ha than redder ones. It is clear that 
the majority of the galaxies detected in Ha have a ((7 — 
B) < 1.0 and that the fraction of detections increases in 
bluer galaxies. In the bluest bin there are more galaxies with 
detected Ha than without, but for {U — B) > 1.0, there are 
more non-detections than detections. 

Previously, Conselice et al. (2007, C07) investigated 
where on the {U — B) verses Mb colour- magnitude di- 
agram galaxies in the POWIR survey fall as a function 
of mass. They find evidence that galaxies evolve onto the 
red-sequence at 2 < 1, but that galaxies with M, > 
10^^'^ M© are generally always red, with a ~40 per cent frac- 
tion that are blue at redshifts greater than z ~ 1.3. 

Figure [5] shows galaxies in this work over-laid on the 
colour-magnitude diagram for the entire POWIR survey, di- 
vided into those with detected Ha emission, and those with 
Ha upper limits (no detection). Of the eight galaxies that 
were targeted and have continuum detections that fall within 
the red sequence, only one has detected star formation in 
Ha. The majority of galaxies without Ha detections lie in 
the red cloud or green valley. Dust will move galaxies to 
redder colours as it attenuates the blue end of the spectrum 
more and so will complicate this plot. 



3.2 Stacked spectra 

To further investigate the galaxies where Ha emission is not 
detected, the spectra from the galaxies without detected Ha 
are stacked (Figure There is no detected Ha emission 
in this averaged spectrum above the noise which equates 
to a mean SFR of 6.6 M0 yr~^ (not dust corrected). This 
value is lower than some individual galaxies with detected 
Ha emission. The non-detections have a mean redshift of 
0.85 while the detections with a SFR less than 6.6 all have 
a redshift less than 0.81. The mean redshift of the galaxies 
with detected Ha emission is 0.78 and this difference in mean 
redshift is partially why the mean SFR limit for the non- 
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Figure 6. a) Fraction of galaxies detected in Ha as a function of redsfiift. There is no obvious bias against detecting Ha emission at high 
redshift in our study, b) Fraction of all our galaxies detected in Ha as a function of stellar mass. Galaxies with stellar masses greater 
than M« = 10^^ Mq have a significantly lower fraction of Ha detections than lower masses, c) Fraction of galaxies detected in Ho as a 
function of U-B colour. There is a significant drop in Ha detections towards redder galaxies. The points have been shifted by ±0.005, 
±0.01 and ±0.005 respectively on the x-axis of each graph for clarity. 
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Figure 8. {U — B) — Ms colour-magnitude diagrams for galaxies detected in Ha (a), and for those not detected (b). The gray points 
are the entire POWIR survey, a) Detections — triangle points have SFR//^ uncorr < SMq yr~^, diamond jjoints have SFR/f„ uncorr 
of 5-15 Mq yr~^ and square points have SFKHa.uncorr > 15 M© yr~^. b) Non Detections - triangle points have SFRHa,uncorr upper 
limit < 5M0 yr~^ and diamond points have SFRjj Q,_„„(,orr upper limit of 5-15 Mq yr~^ 
Ha detections in the red sequence of galaxies. 



SFRs are not corrected for dust. There are no 
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Figure 10. The averaged luminosity spectra for the lowest stellar 
mass bin, M, < ^ q£ q^j. sample. 
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Figure 11. The averaged luminosity spectra for the bluest two 
colour bins in {U — B) for the whole sample, {U — B) < 0.8 (top) 
and 0.8 ^ {U ~ B) < 1.0 (bottom). 



detections is higher than the SFR in some of the galaxies 
with detected Ha emission. 

Figure [To] shows the stacked spectra of the lowest stel- 
lar mass bin when the whole sample (detections and non- 
detections, not dust corrected) is stacked in stellar mass 
bins. There is no significant change in Ha emission in dif- 
ferent stellar mass bins; these emission lines correspond to 
mean SFRs of 9.4 ± 3.3 Mq yr"\ 8.0 ± 3.9 Mq yr"\ and 



Figure 12. Comparison of the dust corrected Ho-derived SFRs 
with the dust corrected UV-derived SFRs for galaxies with UV 
measurements. The triangular points with arrows are upper limits 
for the SFRHa,corr- There is good agreement between the two 
different measures of SFR. 



12.1 ± 5.3 M0 yr"^ in the M. < 10^°'^'^ Mq , lO^^-^^ 
M./Mq < 10" ° and M, > 10"-° Mq mass bins. 

Figure [11] shows the stacked spectra (bot dust cor- 
rected) for the bluest two bins when all the galaxies stacked 
in ([/ — B) colour bins. There is a clear drop in emission 
towards the redder bins corresponding to mean SFRs of 
19.9 ± 3.5 Mq yr-\ 14.0 ± 4.2 Mq yr-\ < 7.3 Mq yr"^ and 
< 6.3 Mq yr-i in the [U ~ B) < 0.8, 0.8 <([/-£)< 1.0, 
1.0 ^{U ~ B) < 1.2 and {U - B) > 1.2 colour bins. 



3.3 Dust corrected star formation rates 

Figure [121 shows the comparison between the dust corrected 
Ha-derived and dust corrected UV-derived SFRs for the 
galaxies with UV measurements. There is good agreement 
between the two SFR tracers. Technically, for each galaxy, 
both SFR tracers have been dust corrected with the same 
E{B — V) value, however, due to the dust attenuation curve 
being used, the actual amount both tracers have been ad- 
justed by is different. As there is a good agreement between 
the SFR measures, this helps confirm the accuracy of both 
measurements (see Bauer 2011). 

Figure 113b . shows the dust corrected Ha-derived SFRs 
plotted against the 24 /im derived SFRs for the galaxies with 
available 24- /im measurements. It is expected that the Ha 
luminosity from star formation correlates with the 24- /xm 
luminosity produced as a result of star formation. However, 
other effects which also produce 24-/im emission would lead 
to an over prediction of the 24 /xm-derived SFR, such as an 
AGN, which however we do not find. 

Figure [T5b shows a relatively good correlation between 
the two measures of SFR apart from two points that lie with 
significantly larger SFRfra^corr than SFR24Mm- We quan- 
tify this difference by the parameter T = (SFRna.corr — 
SFR24;jm)/SFR24,jm- It is uot easy to see how points with 
a high (positive) T are physically possible and they are prob- 
ably due to an overly large dust correction. 

It has also been suggested in some high redshift stud- 
ies that the correction to the colour excess for stellar hght 
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Figure 13. a) Comparison of the dust corrected Ho-derived SFRs with the 24 /^m-derived SFRs for galaxies with 24 measurements. 
Where the extinction correction has been derived from the E(B — V)-SFR relation (as opposed to calculating a galaxy's exact correction), 
the points are plotted in grey. Underneath is plotted the SFR^f^^ excess T, showing 2 galaxies with more star formation derived from 
Ha than 24 /.im measurements, b) shows the same but without the E{B — V)ate.llar = OAAE^B — V)gas correction being applied 



(Equation O does not hold at high redshift, as it over- 
predicts Ha-derived SFRs with respect to the UV-derived 
SFRs (e.g. Erb et al. 2006, Hayashi et al. 2009, Forster 
Schreiber et al. 2009). The correction comes from the ob- 
servation that the Balmer line emission is more attenuated 
(by a factor of ~ 2) than starlight (see e.g. Petrosian, Silk & 
Field 1972, Spitzer 1978, McKee & Williams 1997). When 
models are compared to observations, this effect appears to 
be due to the geometrical distribution of the attenuating 
dust in relation to the ionising stars (Calzetti et al. 1994, 
Chariot & Fall 2000). The fact that it may not hold at higher 
redshifts could be due to different dust geometries, or that 
galaxies at z ~ 2 have a more top-heavy initial mass function 
than the IMF used (see for example Hayashi et al. 2009). 

Figure [T3b plots the Ha-derived SFRs when we do not 
apply this correction (i.e. when using E{B — V)stenar ~ 
E{B — V)gas) against the 24-/im derived SFRs. These 
SFR_f/ci, corr values are lower and the correlation is now closer 
to what would be expected if 24-//m emission comes from 
both star formation and other sources. This is similar to 
what is found in Erb et al. (2006) and Hayashi et al. (2009). 
However Hayashi et al. conclude there is no firm evidence 
that E{B - V)steUar = 0.44 £(5 - V)gas is incorrect. It 
is not yet possible to measure the spatial dust content for 
galaxies at high redshift and thus there is no obvious way to 
validate any theory that the dust geometry in z ~ 1 galaxies 
is different than in local galaxies. We therefore use the 0.44 
correction to the E{B — V)ga3 values in Equation [B] for the 
remainder of this paper. 



3.4 AGN heating and IR emission 

Galaxies with an AGN are believed to heat dust close to the 
engine to a few thousand degrees (see for example Alonso- 
Herrero, 2001) which leads to an increase in the 24 /xm IR 
emission. Therefore, dusty AGN would be emitters of 24- 
^m radiation that could mimic highly star-forming galax- 
ies. X-ray emission is used to distinguish AGN from normal 
galaxies in the EGS field due to both the emission from the 
relativistic jets and from accretion onto the central black 
hole. While X-ray emission is a clear indicator of the pres- 



ence of an AGN, the other types of emission can be at- 
tenuated significantly by dust. This could lead to a large 
population of dusty AGN not detected in X-rays (see for 
example Wang & Jiang, 2006). At higher redshift, the lower 
received fluxes from these AGN could lead to more AGN 
being miss-classifled as star-forming galaxies than at lower 
redshifts. 

To tell whether a galaxy contains an obscured AGN, the 
galaxy's 24-/im emission must be compared to a SFR tracer 
that is not strongly affected by the presence of an AGN. 
Ha is relatively insensitive to temperature and density of 
gas (Osterbrock 1989) and so is only slightly affected by the 
presence of an AGN (e.g. Dong et al. 2008) and thus a better 
representation of the SFR of the galaxy. Figure 113b shows 
the comparison of the 24-/im emission and dust corrected 
Ha SFRs and SFR upper limits for the 18 galaxies that 
have an IR measurement, from the EGS fleld. As the 24- 
/im emission can be caused by star formation, old stars or 
the presence of an AGN, it is expected that most galaxies 
with an AGN would have a higher 24 /xm-derived SFR than 
a dust corrected Ha-derived SFR. There are three galaxies 
which have Ha-derived SFR upper limits that lie below the 
one-to-one relation line and could possibly have a 24-/im 
excess. 

There are four galaxies that appear to have SFRi/a 
higher than SFR24;jm, two of these also have surprisingly 
high Ha/H/3-derived E{B — V) values. These large Balmer 
decrement-derived values may be due to a larger than ex- 
pected [Nil] /Ha ratio. By assuming the UV-derived value 
of E{B — V) is correct and comparing it with the Balmer 
decrement-derived value, a [Nil] /Ha ratio that would be re- 
quired to reconcile the results can be obtained. These values 
are shown in Table |3] for all the galaxies with available H/3 
data. 

The standard use of this type of line ratio diagnostic 
is carried out when plotted against the [OIII]/H/3 ratio, 
for example by Kauffmann et al. (2003b) to separate AGN 
from normal galaxies. The two galaxies that have a positive 
log([NII]/Ha) value would almost certainly be AGN using 
this classification; the three other galaxies appear to have 
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ID log.o (^) 

12024309 0.48 ± 0.32 
12024440 -0.61 ± 0.09 
12024445 1.2 ± 1.8 

12024528 -0.90 ± 0.15 
13035323 <-0.50 



Table 3. [NII]/Ha ratio for galaxies with H/3 detection in the 
DEEP2 spectra. Galaxy 13035323 only has an upper limit. 

log([NII]/Ha) < —0.5 and so probably fall into the star- 
forming region of the diagram. 

The Balmer decrement-derived E{B—V) values are sub- 
ject to large uncertainties due to the modest signal-to-noise 
ratio for the DEEP2 measured emission lines at this red- 
shift. The AGN fraction in the universe at these redshifts 
has been measured to be a few per cent (Silverman et al. 
2009; Haggard et al. 2010) and so the number in a sample 
of star forming galaxies that have been solely selected by 
stellar mass will statistically be low. It is still possible that 
they are highly obscured AGN but there is no evidence from 
our data that this may be the case, and thus we conclude 
that the AGN contamination of our sample is very small. 

3.5 Star formation rate relations 

Figure [141 shows the uncorrected and corrected Ha SFRs as 
a function of stellar mass. In both cases there is a possible 
slight increase of star formation with mass but is otherwise 
constant at all masses at M, > 10^*^ M© . The average points 
only use the measured SFRs, and do not include the upper 
limits, and these averages stay roughly constant with stellar 
mass. 

Previous observations have suggested that there is a 
"Galactic Main Sequence" (GMS) of star forming galaxies 
on a SFR - M, diagram (Noeske et al. 2007a hereafter N07, 
Daddi et al. 2007a hereafter DQ7, Pannella et al. 2009) and 
two of these lines are overplotted. The N07 line (dots) was 
measured for a sample of galaxies in the redshift range 0.2 < 
2 < 0.7, and has a slope of log (SFR) oc (0.67 ± 0.08) log M. . 
They suggest that at redshift one, their measured GMS 
would appear to have the same slope but lie ~ 3 times higher 
than their line for 0.2 < 2; < 0.7. Their GMS appears to 
lie below the majority of our galaxies, suggesting that our 
galaxies have higher measured SFR at a given stellar mass. 

The line from D07 (dashed) has a slope of log (SFR) oc 
0.9 log M« which is similar to Forster Schreiber et al (2009) 
who found it to be close to their Ha derived values. The 
line lies close to our dust-corrected detections, however it 
appears to have a steeper slope than our galaxies. 

We use a Monte Carlo simulation to determine the prob- 
ability that our data fits the N07 GMS. Each data point is 
moved randomly within its error bar and a line of best fit 
is found using a least-squares method. This is repeated 10* 
times and the mean slope and intercept point are calculated 
as: 

\og{SFRHc^corr) = (0.40 ± 0. 13) log M. - (2.4 ± 1 .3) , (10) 
suggesting a slope slightly flatter than found in previous 



work. There is a 8 per cent probability that the slope for 
our data is within the errors of the slope of N07. There is no 
confldence greater than the 1 cr level that our derived slope 
is different to the slope of N07, so we cannot confirm that 
the two are the same. 

A measure of the formation histories of galaxies can be 
inferred from the specific star formation rate (SFR divided 
by stellar mass, SSFR) which is plotted against mass in Fig- 
ure llSI The SSFR takes into account the fact that the same 
SFR will have a different effect on galaxies with different 
stellar masses and so gives a clearer picture of the SFR his- 
tory of the galaxies probed. 

The SSFR appears to fall with increasing mass, espe- 
cially at M« > 10^"'^ M0 , and the mean value of the detec- 
tions is possibly constant at masses lower than this. Along 
with the previous finding that there are no galaxies with 
M» > 10^^ Mq that have Ha detected at a redshift below 
z = 0.93 in the sample, this suggests that at z ~ 1, galaxies 
with M* > 10^^ Mq are much less likely to still be forming 
stars; but where they are, they are forming stars at a similar 
relative rate to lower mass galaxies. 

Overplotted on Figure 115b is the theoretical relation 
from the exponential star formation history model of Noeske 
et al. (2007b) at a redshift of z = 0.8. This model uses a 
mass dependent formation redshift for galaxies and evolves 
their star formation using a closed-box model. The model 
appears to fit the relation observed in our data. However, 
the dust corrected SSFRs for galaxies with Ha detections 
are possibly slightly higher. When taking into consideration 
the upper limit points, the model is a slightly better fit. Our 
points do though lie within the scatter of the data used to 
derive this relation in Noeske et al. 

In Figure [161 the SFR is plotted against {U — B) colour 
along with the mean value of the detections. Also plotted 
is the mean colour of all the points. As the number of Ha 
detections drop with increasing [U — B) colour (Figure |6t), 
the number of galaxies with Ha detections is very low at 
{U — B) > 0.9. The galaxies that are detected are the galax- 
ies with the most on-going star formation and so are not 
necessarily representative of all galaxies. Taking the mean 
of all the points (detections and upper limits) does not nec- 
essarily give a meaningful number but it better shows the 
overall trend alluded to by the data points. It also allows 
the large 0.9 ^ {U — B) ^ 1.4 bin to be split into two. It is 
clear that the two methods show different trends with the 
mean value of detections constant and the mean value of all 
points dropping for redder galaxies. 

Comparing the two methods, there are apparently two 
populations of galaxies, those that are star forming which 
have roughly the same SFR whatever their colour; and 
galaxies that are not star forming. At redder colours, there 
are more galaxies not star forming which lowers the average 
SFR of all the galaxies at very red {U ~ B). 

To look at these effects more productively, the SSFR is 
compared to the {U - B) colour in Figure [TTl The SSFR is 
clearly lower in redder galaxies and the trend is similar for 
both the mean value of detections, and the mean value of 
all points. 
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Figure 14. The uncorrected (a) and corrected (b) Ha SFR plotted against stellar mass for all galaxies. Galaxies with Ha detections are 
plotted as crosses, the galaxies with an Ho upper limit are plotted as a triangle with a downward pointing arrow. Where the extinction 
correction has been derived from the E{B — y)-SFR relation (as opposed to calculating a galaxies exact correction), the points are 
plotted in grey (plot (b) only). The large circled stars are the mean of the detections in mass bins and the best fit line for the galaxies 
with detections is plotted as the dashed line on each plot. The positions of the galactic main sequence as measured by Noeske et al. 
2007a (NOT, dotted) and Daddi et al. 2007 (D07, dot-dashed) are shown, as well as the best fit line for our galaxies with Ha detections 
(dashed). 
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Figure 15. The uncorrected (a) and corrected (b) Ha SSFR plotted against stellar mass for all galaxies. The points are the same as in 
Figure [Ml The solid line is the modeled relation for galaxies with an exponential star formation history, as discussed in Noeske et al. 
(2007b). 




0.6 0.8 1.0 1.2 1.4 0.6 0.8 1.0 1.2 1.4 

U-B U-B 



Figure 16. The uncorrected (a) and corrected (b) Ha SFR plotted against colour for all galaxies. The points are the same as in Figure 
1141 The large circled stars are the mean of the Ha detections in mass bins, the square points are the mean SSFR of all galaxies in mass 
bins, including the upper limits — See text for details. The circled stars and square points are shifted by ±0.005 in the x axis for clarity. 
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Figure 17. The uncorrected (a) and corrected (b) Ha SSFR plotted against {U — B) colour for all galaxies. The points are the same as 
in Figure [T6l 




Figure 18. a) CAS asymmetry against dumpiness parameters for the galaxies with 24-^m measurements. The parameter space for 
mergers is defined as A > 0.35 and A > S. The points are sized by the SFR excess T. b) Concentration against asymmetry for the 
same galaxies. Mergers are defined as A > 0.35. The boundary between Early and Mid/Late type galaxies is C = 21.51og^Q(A) + 31.2 
(ConseUcc et al. 2008b). 



3.6 CAS morphology 

The CAS morphology parameters (Conselice 2003) are anal- 
ysed on our sample using the ACS imaging to determine 
whether any of our sample are likely mergers which may af- 
fect the SFR determinations. These parameters also allow 
us to examine the nature of the star formation we measure 
and investigate the morphologies of the galaxies. 

Within the asymmetry (A) against dumpiness (S) plot, 
a galaxy merger is defined as having A > 0.35 and S > 
A (Conselice, Rajgor & Myers 2008b). For a concentration 
(C) against asymmetry plot, a galaxy merger is defined as 
having A > 0.35 and C < 2.44 log( A) -1-5.49. (Conselice et al. 
2008b). These plots are shown in Figure [181 along with the 
position of the three galaxies identified from the imaging; 
the points are coloured by bins of SFRi/a excess T. 

Plotted on Figure 118b is the line dividing early- and 
mid/late- type galaxies as defined in Conselice et al. (2008b). 
It is clear that all of the galaxies we observe that are not 
probable mergers, are mid- or late-type spirals which is con- 
firmed by the visual morphology (Figure [21). 

There is no clear relation between the SFR/^a excess 
T and the position in the CAS parameter space. The SFR 



excesses are therefore likely due to the individual circum- 
stances of each galaxy. 

The possible merging galaxy (12024436) is confirmed to 
be a merger by the CAS parameters. Two other galaxies lie 
very close to the boundary of the merger parameter space 
- galaxies 12024445 where the upper error bar extend into 
the merger region, and 13035323. Galaxy 13035323 does not 
have an unusual SFRjfa.corr however galaxy 12024445 is one 
of the galaxies with a larger SFR_f/a,corr than SFR24;jm. 

The possible merging systems 12020031 and 13026831 
are not mergers as found in CAS parameter space. The dou- 
ble galaxy system 12020027 does not appear to be a merging 
system likely because the CAS code ignored the other galaxy 
(see Conselice 2006 for time-scale issues within CAS). 

Figure [191 shows the asymmetry and dumpiness of our 
galaxies as a function of their measured SFRs. There is pos- 
sibly a relation between the asymmetry and the SFR and 
the galaxies are fit by a line with the relation 

A = (1.09 ± 0.36) X 10"^ SFRnc^.corr + (0.182 ± 0.034) (11) 

This is likely to be due to the star formation happening 
in an asymmetric manner - the more star formation hap- 
pening, the more asymmetric the galaxy is - or that what 



16 



Twite et al. 




200 
(M. yr-') 



300 




200 
(M. yr-') 



300 



Figure 19. a) CAS Asymmetry against corrected SFR for the same galaxies as Figure lTsl witii tlie best fit line sliown. b) CAS dumpiness 
against corrected SFR for the same galaxies. 



ever mechanism is causing the asymmetry is driving the star 
formation. If this were the case, it would point to star for- 
mation not being 'triggered' by some (possibly small) event 
that then causes a large 'snowballing' of star formation ac- 
tivity, but a regime where larger events lead to larger star 
formation bursts. In reality both effects probably happen ac- 
centuating the effect: star-formation triggering events cause 
the galaxy to look more asymmetric and trigger asymmetric 
star-formation episodes. 

Somewhat surprisingly then there is no obvious relation 
between a galaxy's SFR and its dumpiness (Figure 119b ). 
This might imply that star formation is happening in large, 
less well defined regions of the galaxy. There is no relation 
observed between the concentration and SFR. 



4 DISCUSSION 

4.1 The shut-ofF of star formation in galaxies 

In the scenario of downsizing the most massive galaxies 
should turn off their star formation at the redshift ranges 
we probe. In our sample, the most massive galaxies with 
star formation and with stellar masses M. > 10" M0 are 
found only at a redshift greater than z — 0.93. The drop-off 
in the number of star forming galaxies with masses above 
M* — 10^^ Mq is sharp, suggesting that the process of turn- 
ing off star formation happens very quickly, in less than a 
Gyr. This appears to be a sudden quenching rather than a 
slow ending of star formation, similarly it is not likely to be 
a "running out" of gas within star forming regions as in this 
scenario we would expect this to lead to a declining SFR 
over a longer period of time. 

There is only one galaxy that is detected in Hq and lies 
within the red sequence of galaxies in our sample, show- 
ing that star forming galaxies detected with Ha at z~l 
are mostly blue. There are only a small number of non- 
detections (4 out of 25) that lie definitively within the region 
of blue galaxies and these could be either non-detections of 
star forming galaxies or actual non-star forming galaxies. 

In line with what other studies have found at other 
wavelengths, this study finds a slightly increasing Ha- 
derived SFR with stellar mass, and a decreasing SSFR with 
mass relations. There is some suggestion that the normal- 



isation of the 'galactic main sequence' (GMS, N07) is de- 
pendent on the wavelength and method used to measure 
the SFR, and the Ha-derived SFRs appear to be slightly 
higher than other measurements. This could be the result 
of the sensitivity of measuring SFRs which in this study are 
limited to greater than a few Mq yr~^. With a dust correc- 
tion, this becomes a lower limit of SFRs greater than a few 
tens of Mq yr~^. [Oil] star formation rates can be measured 
(uncorrected) typically down to 0.5-1 M0 yr~^, which when 
corrected give SFRs with typically half the value detected in 
this sample. However, to reconcile our results with the N07 
derived GMS, there would have to be a significant number 
of galaxies which have a SFR an order of magnitude below 
the level of our GMS at all masses. 

There are two probable causes of our detections hav- 
ing larger average SFRs than the GMS measured by N07. 
A major contributor will be that some of our galaxies have 
been selected to have significant 24-/im emission, meaning 
that they are either highly star forming, or have a signifi- 
cant AGN component. Two galaxies (12024309 & 12024445) 
could have [Nil] /Ha ratios indicative of AGN, but have a 
SFRhq excess suggesting otherwise. These two therefore 
probably have a slightly larger amount of [Nil] emission 
leading to a overestimation of the Ha-derived SFR. There 
is little other evidence of any other galaxies in this sample 
having a dominating AGN and so the higher SFRs mea- 
sured are likely because the galaxies actually have a higher 
rate of star formation. These would then lie above (i.e. with 
a higher SFR) any GMS relationship found in more com- 
plete samples. The galaxies however lie much above the N07 
GMS (even with their estimated raising of the normalisa- 
tion for higher redshifts), significantly beyond their 0.3 dex 
estimated scatter. 

The Daddi et al. (2007, D07) GMS is derived from UV 
and 24-/im emission(as opposed to [Oil] and 24-/im in N07) 
and appears to be at a similar level to our data, but it is not 
clear whether our data follows their slope. It is clear then 
that the choice of SFR indicator can be important. C07 mea- 
sured that the 24 /xm-derived SFR were considerably higher 
than the [Oil] values and it is not unusual that the two 
methods do not give agreement on a galaxy-by-galaxy basis 
(e.g. Salim et al. 2007). 

The SSFR is a measure of how much of a galaxy's mass 
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ID E{B - V) SFRhccott T C a S Notes 



12020027 


1.02 


101 ± 31 


10.6 


2.70 ± 0.09 


0.15 ± 0.02 


0.02 ± 0.002 


Two galaxies in image 




12020031 


0.68 


44 ± 23 


1.84 


2.96 ± 0.03 


0.33 ± 0.05 


0.47 ± 0.02 


Possible second galaxy in 


image. 


12024436 


0.94 


124 ± 55 


0.472 


2.65 ± 0.03 


0.66 ± 0.03 


0.13 ± 0.01 


Possible star in slit. CAS 


merger. 


12024440 


0.60 


11.1 ± 6.4 


1.63 


3.12 ± 0.07 


0.16 ± 0.03 


0.06 ± 0.003 






12024445 


0.85 


136 ± 30 


4.20 


3.15 ± 0.04 


0.35 ± 0.05 


0.04 ± 0.004 


Possible CAS merger. 




13026831 


0.66 


52 ± 37 


1.10 


3.81 ± 0.17 


0.20 ± 0.01 


0.10 ± 0.002 


Poss merger in image, not 


in CAS. 


13035323 


0.64 


32 ± 19 


0.306 


3.38 ± 0.07 


0.33 ± 0.02 


0.10 ± 0.004 


Possible CAS merger. 





Table 4. Galaxies with high Ha-derived SFR excess, T, have morphological features, or identified as mergers by CAS parameters. 



has been formed in recent star formation episodes. The most 
massive galaxies with M* > 10^^ Mq have a lower SSFR 
than the less massive galaxies which have a roughly constant 
value. These very massive galaxies have had less of their 
mass formed in recent star formation and so must contain 
older stellar populations. This is a clear indication that the 
most massive galaxies had their star formation quenched at 
a higher redshift than lower mass galaxies. 

4.2 Colour transition 

We find a separation between the star forming and non-star 
forming galaxies in SFR-colour plots (Figure [16]). Gala^des 
that are towards the red end of our sample (in {U — B)) 
and have measured star formation, have the same average 
SFR as the blue galaxies. In the SSFR-colour plots (Figure 
I17p . the SSFR drops with colour due to the fact that the 
red star forming galaxies are more massive than the blue 
galaxies with the same SFR. This is seen in Figure [20] and 
is possibly a selection effect. Half the galaxies in the sample 
were selected to have a 24- /xm detection - for a galaxy to 
have enough star formation so that it gives 24-/im emission, 
but not be significant enough to make the galaxy blue, the 
galaxy must be more massive. Figure [20] also shows a consid- 
erable difference in the slope of this relation for star-forming 
and passive galaxies. The star-forming galaxies are fit by a 
line with the relation 

log (M.) = (1.87 ± 0.33) X {U -B) + (9.23 ± 0.28), (12) 
and the galaxies without Ha detections are fit by the relation 

log (M,) = (0.61 ± 0.27) X [U -B) + (10.27 ± 0.28). (13) 

A two dimensional Kolmogorov-Smirnov test (K-S test; 
Fasano & Franceschini, 1987) indicates that these two pop- 
ulations are however not significantly different. 

A more detailed look at the SFR and SSFR against 
colour relations can be obtained by binning the points by 
stellar mass (Figure [2T|) . The effects of downsizing are clearly 
seen in these plots with the highest mass bins having many 
more passive galaxies than star-forming ones, similarly the 
lowest mass bin has more star-forming than passive galax- 
ies. Both the lowest and intermediate mass bins in the SFR- 
colour plot can be described as a number of star-forming 
galaxies with a range of SFRs but with similar colours, with 
the passive galaxies (with similar SFRs) having a range of 
{U — B) colours. The most massive galaxies are mostly pas- 
sive, red galaxies. 




(U-B) 

Figure 20. Stellar mass against {U — B) colour for galaxies di- 
vided into those with detected star formation (crosses) and no 
detected star formation (triangles). The redder galaxies tend to 
be more massive than the blue galaxies — this is possibly a se- 
lection effect due to half our galaxies needing 24 /xm emission for 
selection and thus have star formation. There is a lack of massive, 
blue galaxies which should have been detected. The solid line is 
a least-squares fit to the galaxies with Ha detections and has a 
slope of log(M,) oc (1.87 ±0.33) x (U - B). The dashed line is 
a least-squares fit to the galaxies without detected Ho and has a 
smaller slope: log (M*) oc (0.91 ± 0.18) X (U - B) 

When galaxies turn off their star formation, they also 
usually start changing colour from blue to red. This would 
lead to the galaxies moving downwards and to the right in 
the SFR-colour diagram. If the concept of downsizing is cor- 
rect, then the most massive galaxies shut off their star for- 
mation before less massive systems. We can then use a mass 
cut to separate a population of galaxies that have shut off 
their star formation (the most massive galaxies) from a pop- 
ulation that are yet to shut off their star formation (the least 
massive galaxies) . Figure [21] splits the galaxies into stellar 
mass bins and shows their position on a SFR-colour diagram. 

In each stellar mass bin, the star forming galaxies have a 
bluer {U — B) colour than the passive galaxies and the range 
of colours that the star forming galaxies have in each stellar 
mass bin is quite narrow. Passive galaxies have a range of 
{U — B) colours. Each stellar mass bin is successively redder, 
this is due to the relation between stellar mass and colour 
(Figure !^ . 

As seen in Figure [T7] there is a relationship between the 
SSFR and colour. This is partially expected - SSFR traces 
the effect that star formation has had on the the recent his- 
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Figure 21. SFR vs colour binned by stellar mass. Filled in symbols have detected Ha emission, empty symbols do not. The black points 
are the galaxies within the stellar mass range of the bin: Triangles - M, ^ 10^0 '' M© , diamonds - 10^0 '' Mq M* ^ lO^^ O Mq , 
squares — M« ^ lO^'^" Mq ; the light grey points are the galaxies in the other mass bins shown for comparison. The lowest mass bin 
contains the highest fraction of star forming galaxies, the highest mass bin contains the highest fraction of passive galaxies. Within each 
mass bin, galaxies that are star forming have similar (U — B) colours and a range of SFRs; galaxies that have finished their star formation 
have a wide range of colours. 



tory of the galaxy. This will also be borne out in the colour 
of the galaxy, with the more recent the star formation, the 
more affect it has on a galaxy's colour. This then explains 
why the different mass bins are separated along the colour 
axis of the SFR-colour diagram (Figure I^Tj) : the redder sys- 
tems have lower SSFRs and so for the same SFR they must 
be more massive than bluer systems. 

This relation could be affected by our selection criteria. 
We could be missing low-mass passive objects which are red. 
However due to their low mass, they would not lie consid- 
erably below galaxies already plotted on the SSFR-colour 
diagram. 

The galaxies investigated appear to follow what would 
be expected in a downsizing senario. There are significantly 
fewer galaxies with M* > 10^^ Mq forming stars and the 
most massive galaxies are only forming stars at the highest 
redshifts probed. More detailed studies into the SFR-M* as 
a function of redshift are not possible due to the small size 
of our sample. 



5 CONCLUSIONS 

We have analysed NfR spectroscopy of 4f massive (M« > 
10^° Mq ) galaxies taken from the POWIR Survey of which 
half were chosen to have significant 24-/im flux, no X-ray 
emission and within the range 0.4 < Zspec < 1-3. We have 
acquired NfR spectroscopic data for these systems and mea- 
sured SFRs from the Ha emission line. The Ha emission line 
is the best direct measure of instantaneous star formation 
as it is coupled directly to the massive stars in formation, is 
largely independent of metalicity and is less effected by dust 
attenuation than either the [Off] emission line or UV flux. 

We investigated several methods of correcting high- 
redshift galaxies for dust using multi-wavelength methods 
and chose to ultimately apply a method derived from UV 
data. We flnd that the Hq SFR measurements are in good 
agreement with those measured by UV emission, but are 
slightly higher than those measured from emission. 
By not using the correction to the E{B — V) value due 
to differences associated with gas emission and stellar hght 
(Calzetti et al. 2000) our Ha corrected SFRs are similar to or 
less than the galaxies' 24-^m derived SFRs. The AGN con- 
tamination of our sample is low - although highly obscured 



AGN may exist within the sample there is no evidence that 
there are any. 

We find that there is a drop in the fraction of galaxies 
with detections of Ha for systems with a stellar mass greater 
than M. = f O" Mq . The massive galaxies with detected Ha 
are all found at redshifts greater than z = 0.93. The fraction 
of galaxies with detected Ha also drops for galaxies that are 
red in U-B colour. The sharpness of this drop implies that 
star formation is truncated over a very short period of time 
rather than dropping steadily over an extended epoch. 

We investigated the Ha SFR as a function of stellar 
mass and find that SFR stays roughly constant with stellar 
mass. Likewise, the SSFR drops with increasing mass, espe- 
cially at M* > fO^°'^ Mq . This implies that at z ~ f , galax- 
ies with M* > 10^^ Mq are less likely to be forming stars, 
and when they are, they are doing so at a much lower rate 
relative to their mass. These massive galaxies have a lower 
SSFR than lower mass galaxies, implying that they have on 
average a stellar population that formed earlier, and must 
have had their star formation quenched at a higher redshift 
than the lower mass galaxies. 

We find that the observed slope in the SFR-mass rela- 
tion appears to be slightly flatter than what is found in pre- 
vious studies. The SSFR-mass relation also agrees with the 
relation predicted by an exponentially decreasing star forma- 
tion history model with mass dependent formation redshifts. 
Our measured SSFRs however are slightly higher than the 
Noeske et al. (2007b) model predictions. 

We also investigate the SFR as a function of colour and 
flnd that the galaxies that are detected in Ha form stars at 
similar rates independent of colour. When split into bins of 
stellar mass, within each bin we flnd star forming galaxies 
that have a narrow range of {U — B) colours while galax- 
ies not detected in Ha have a wider range of colours which 
are mostly redder than the star forming ones. These mass 
bins can be viewed as populations of galaxies in different 
stages of their evolution from star forming to passive, and 
suggests that galaxies drop in star formation before becom- 
ing redder. We also flnd that the speciflc SFR drops with 
increasing {U — B) colour (redness) for all galaxies; this is 
probably due to both being a measure of the effect on the 
galaxy of recent star formation. This relation leads to the 
more massive galaxies (with lower SSFR) becoming redder. 
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Evolutionary tracks on the SFR-colour plot are therefore 
shifted towards redder colours for the most massive galax- 
ies. We also investigate the structures of our sample which 
have ACS HST imaging. We find that there is possible re- 
lation found between a galaxy's SFR and its asymmetry as 
measured by the CAS parameters (Conselice 2003). 

Ha-derived SFRs at 2: ~ 1 appear to be similar to those 
measured using other tracers. There is a large body of work 
studying star formation at these redshifts using these other 
tracers, and we have shown, using the more reliable Ha val- 
ues that these other results are consistent at least up to 
z ~ 1.4. As NIR spectroscopic instrument technology con- 
tinues to evolve, especially on 8-lOm telescopes, the growth 
in the number and quality of Ha measurements, particularly 
at high redshift, will allow us to compare these different SFR 
tracers at higher redshifts and therefore obtain an even more 
accurate measure of the star formation history of galaxies. 
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